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Abstract 

A time-evolution of quantum meson fields is investigated in a linear sigma model 
by means of the time- dependent variational approach with a squeezed state. The 
chiral condensate, which is a mean field of the quantum meson fields, and quantum 
fluctuations around it are treated self-consistently in this approach. The attention 
is payed to the description of the relaxation process of the chiral condensate, where 
the energy stored in the mean field configuration flows to the fluctuation modes. It is 
shown that the quantum fluctuations plays an important role to describe this relaxation 
process. 
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1. Introduction 



One of the recent interests associated with the physics of the relativistic heavy-ion col- 
hsions is the investigation of the dynamics of matter at very high energy density. As for 
the chiral phase transition, especially, the possible formation of a disoriented chiral conden- 
sate in ultra-relativistic nuclear collisions has been discussed. "0) The time-evolution of 
the chiral condensate in the chiral phase transition is a main subject to study in diverse 
work, namely, the dynamics of the phase transition is investigated. In these studies, the 
0(4) sigma model has been used mainly as an effective model of quantum chromodynamics 
(QCD). The time-evolution of a scalar field is also important to describe some aspects of 
inflationary model of universe. "ll^l However, most analytical results in this field concern 
the classical equations of motion. 

When the high energy nucleus-nucleus collisions occur, one expects the formation of the 
quark-gluon plasma. Also, it is expected that the chiral phase transition occurs and the chiral 
symmetric phase is realized. After the collisions, the hot region created by the collision will 
expand spatially and the temperature will become lower. Thus, the chiral phase transition 
will occurs again and the chiral broken phase will be realized. Here, let us consider so-called 
quench scenario. The chiral condensate which is a chiral order parameter must grow from 
to (cr). The relaxation process should be realized since the energy stored initially should 
flow elsewhere. Thus, the interesting and important problem is to make it clear how the 
relaxation occurs. The dissipative process in the linear sigma model has been treated in the 
framework of the Caldeira-Leggett theory & in the context of the decay of the disoriented 
chiral condensate. In their work, the mechanism of relaxation has been given within the 
classical level. 

One of the present authors (Y.T.) together with D. Vautherin and T. Matsui has indicated 
that the quantum fluctuation plays an important role for the chiral phase transition^ and 
the relaxation of the collective isospin rotation of the chiral condensate.© The Gaussian 
approximation of the functional Schrodinger picture 111' for the field theory has been adopted 
in our previous work. In this approach, both the mean field and the quantum fluctuations 
can be treated self-consistently including the quantum effect of the higher order of h. This 
approach is essentially same as the time-dependent variational approach with a squeezed 
state. In this paper, we will use the later approach to investigate the time-evolution of 
the chiral condensate and the quantum fluctuation around it because the fluctuation modes 
can be constructed explicitly. Except for the previous work,0'''Ei'' the time-development 
of the chiral condensate was examined to leading order in a large expansion in Ref.^ 
by Cooper et. al. However, in our time-dependent variational approach with a squeezed 
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state, it is easily possible to investigate the role of the quantum fluctuations systematically 
because the quantum state is definitely constructed by means of the squeezed state and 
the mode-expansion for quantum fluctuations is explicitly performed. The main aim of this 
paper is to point out that the quantum fluctuations are responsible for the occurrence of 
the relaxation process of the chiral condensate or chiral order parameter. Of course, the 
expansion of the system leads to a friction term.i^'0^ However, even if the expansion is not 
taken into account, it will be shown that the behavior like the damped oscillation for the 
chiral condensate appears due to the inclusion of the quantum fluctuations. 

This paper is organized as follows : The time-dependent variational approach with a 
squeezed state is reviewed briefly and the mode functions in this approach are introduced 
in the case of the uniform condensate in the next section. This approach is applied to the 
linear sigma model and the basic equations are obtained in §3. The expanding geometry 
is presented in Appendix. The numerical results and discussion are given in §4. The last 
section is devoted to a summary. 



§2. Time-dependent variational approach with a squeezed state 

In this section we briefly give the basic ingredients in the time-dependent variational 
approach with a squeezed state. S The main advantage of this approach to the dynamical 
problem including both the degrees of freedom of the mean field (order parameter) and the 
quantum fluctuations around it ("mesons" if quantized) is in the fact that both degrees of 
freedom can be treated self-consistently through the nonlinear coupling. Further, we can get 
the canonical equations of motion describing their time- development. 

2.1. Basic equations of motion 

Apart from the sigma model, for simplicity, let us first consider a scalar field theory with 
one component. We start with the following Hamiltonian : 

H = j d^x\^-^{xf + \ {V<P{x)f + U[<P{x)]] . (2-1) 

We take a following squeezed state as a trial state : 

= exp (Sit)) ■ Nit) exp (T(t)) |0) , (2-2) 
S{t) =i J d^x[7r{x,t)(f){x) -Tp{x,t)7r{x)] , 

T{t) = J Jd'xd'y 0(x) \^-^-{G-\x,y,t)-G^'\x,y)-')+zE{x,y,t)^ 0(2/), 



3 



where |0) is the reference vacuum with respect to free boson annihilation operators a^; 
afc|0) = 0, and G^'^\x,y) = {0\(p{x)(p{y)\0) . Further, is a normalization factor. The vari- 
ational functions are Tp{x, t), w{x, t), G{x, y, t) and S{x, y, t), respectively. The expectation 
values for the field operators are obtained as 



{<P{tMx)c^{y) \<P{t)) = ^{x, tMy, t) + G{x, y, t) 

Ut)+A{x\mG{t)E{t)\y) . 

(2-3) 



{m\A^My)\m) =n^,tMy,t) + lG-'{x,y,t) + ^x\E{t)G{t)^^^^ , 



where we used the following notation : 

{x\S{t)Git)S{t)\y) =J J d'x'd^y'U{x,x',t)G{x',y',t)U{y',y,t) . (24) 

Here, ^(cc, t) represents the condensate or the mean field. The two-point function G{x, y, t) 
means the two-point correlation function with the order of %. This function represents the 
quantum fluctuation around the mean field Tp{x,t). For later convenience, we introduce 
M^'^^[f{x,t)] defined as 



M^^\Tp{x, t)] = exp ^, t)^} 



(2-5) 

z=Tp{^,t) 



with f/*^")[2;] = d"-U / dz''\ Here, A^*^°^[^(a;, t)] is nothing but the expectation value of the 
functional t/[0(a;)] of the field operator 0(a;). From this expression, it is seen that the 
fluctuation around the mean field can be taken into account with the loop contributions 
G(cc, X, t). 

The state \^{t)) is determined as an functional oiTp{x^ t), W{x, t), G{x, y, t) and U{x, y, t), 
and their time- development is governed by the time-dependent variational principle : 

s J dt{<pm^ - H\m) = . (2-6) 

The expectation value of the Hamiltonian is expressed as 

H^m\H\m)- Jd'xEix), (2-7) 
£{x) = lwix,tf + i iVTpix,t)f + M^'\Tp{x,t)] 

+l{x\G-\t)\x) +2{x\mG{t)m\^) + ^ lim VxVy(x|G(t)|2/) . 
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The variational equation ( p-6| ) gives the canonical equations of motion : 

dip{x,t) _ 6H dw{x,t) _ 5H 

dt 5w{x,t) ' dt 5Tp(x,t) ' 

dG{x,y,t) 6H dE{x,y,t) 6H 



(2- 



dt 6E{x,y,t) ' dt 6G{x,y,t) ' 

As a result, we obtain the following equations of motion : 

- VM^, t) + M'^'^ [Tp{x, t)] = , (2-9) 

{-Vl+M^^\Tp{x,t)]]G{x,y,t) = ^{x\G-\t)\y) -2{x\ij{t)G{t)\y) 

-A{x\E{t)'G{t)\y) , (2-10) 
= 2[{x\G{t)E{t)\y) + {x\E{t)G{t)\y)] , (2.11) 

where we used the equation dt'^ = vr. The above-derived equations of motion are the basic 
equations in the squeezed state approach. 

2.2. Mode expansion 

From now on, let us concentrate our attention on the equations of motion for quantum 
fluctuations in Eqs.(2-10) and (2-11). We assume that the condensate is uniform spatially 
with a translational invariance. With this assumption, the two-point functions G{x, y, t) and 
U{x, y, t) only depend on x — y. Thus it is possible to carry out the Fourier transformation 
as follows : 



G{x, y,t) = j^J d'ke^'^-^^-y^G^it) , y,t) = j^^l ci3fee*-(^-y)rk(t) . 

(2-12) 

From Eq.(2-ll) the equation Gk = 4Gk^k is obtained. By using the above transformation 
and Eq.(2-10), we can get the following equation of motion for the Fourier modes G'k(^) : 

¥-<f^ - Wf) + - 4CW, + 1^'"'!^ = " • f^-"' 

where [TW-^^GJk means the Fourier-transformed mode of Ai^'^^ [Tp{t)]G{x,y,t) : 
M^^^[ip{t)]G{x,y,t) = (l/(27r)3)/rf3A;e^'^-(^-y)[A^(2)G]k. It should be noted here that Gk 
has the meaning of the Gaussian width in terms of the functional Schrodinger picture, 
that is, Gk is positive definite. We can thus express Gk(t) as a square of the mode function 

GkW = vUt? ■ (2-14) 
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Substituting Eq.( |2-14] ) into Eq. (|2-13| ), the equation of motion for fluctuation mode is obtained 



as 

^.W + feVW + IA^'^'Gh^-j^^^O. (2.15) 

The energy density in Eq.(2-7) is expressed in terms of the mean field ^(t) and the mode 
functions as follows : 

Here, we give two comments : One is about the mode expansion. If we rewrite the mode 



function into 77k(i) = 1/ y'^^kit), the two-point correlation function G{x, y, t) is expressed as 
G{x,y,t) = Ek^k{x,t)^l^{y,t)/{2e]^{t)) where ^k{x,t) = e*-''e-^/'^k(*')dt'. rj^^^jg expression 
is nothing but that of the usual mode-expansion by the plane wave. Here ek(t) corresponds 
to the time-dependent single particle energy. Another is the treatment of the system at 
finite temperature. If it is necessary to deal with the system at finite temperature in the 
thermal equilibrium state, the mode expansion must be slightly modified. In the thermal 
equilibrium, the variables are time- independent. The field operator is divided into = Tp + (f). 
Further, is expanded by the mode functions : 

= E^(^k&k + a&l) , (2-17) 
k V^ek 

where 6k I ^) = at zero temperature. In the thermal equilibrium state, we assume that 
(6k&k') = (^k^k') = s-iid (6[.6k') = '"'k'^kk' where rik is the bose distribution function : 
nk = l/(e'^'^k — 1). Then, the two-point function G{x,y) is expressed as 

Gi^^y) = ^^M^my) + m^omiy) + Cki^)^M)nk] . (2-18) 



Thus, it is necessary to replace the Fourier transformation in Eq.( p-14D into 

Gk = {l + 2n^)vl (2-19) 
when we treat the system at finite temperature in the thermal equilibrium state. 

§3. Squeezed state approach to the hnear sigma model 

In this section, let us derive the equations of motion for the linear sigma model in this 
approach developed in §2. The Hamiltonian density is given as 

H = l-TTaixf + IVM^)^M^) + A (M^f -V^Y - hMx) (3-1) 
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where a runs ~ A^ — 1. The index means the sigma- field and 1 ~ A^ — 1 mean the pi- fields. 
Hereafter, we introduce = 4Af ^. The trial state we here adopt is expressed as 



N-l 



= n expiSait)} ■ Kit) exp{r,(t)}|0) . 

a=0 



(3-2) 



The expectation value of Hamiltonian density, (<P(t)|7i|<P(t)), is easily calculated similar to 
Eq.(2-7). For later convenience, we define A^^"^-* same as Eq.( p-5D : 



M^:\if,{x, t)] = exp I Ig,{x, x,t)S^]-^l -^z^ + Xz' 



dz"^ dz^' 



(3-3) 



with a = ~ — 1. We omit the constant term m^/(16A) because of no contribution to 
the dynamics. 



The time-dependent variational principle ( p-6[ ) gives us the canonical equations of motion 
for (Tp^,7Ta) and Ea) and leads us to the following equations : 



dlp,{x,t) _ SH 
dt 

d7Ta{x,t) 



6TTa{x,t) 

SH 



Tra{x,t) 



dt 



^^a{x,t) 

-VVa(^,t)+-Mi'^[^a(^,t)] 

+AXip^{x, t) ^ (^b(a;, tf + Gb{x, x, t) ) -h5, 



(3-4) 



dGa{x,y,t) bE_ 

dt 

dEa{x,y,t) 



dt 



6Ea{x,y,t) 
5H 
5Gaix,y,t) 



2 [{x\Ga{t)Eait)\y) + {x\E^it)Gait)\y)] 



lM^^^[if^{x,tMx -y)- \{^\-^,\y) + '^{A^amv) 

+ \'^x^y5{x -y) + 2X J2 {^bi^, t? + Gb{x, x, t) j 6{x - y) 



(3-5) 



The last equation in Eq.(3-5) is simply expressed by multiplying Ga{x,y,t) from the right- 
hand side : 



{-Vl, + M'^:\ip,{x,t)]}Ga{x,y,t) 
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■(a;|G':i(t)|t/) -2(a;|i:„(t)G'„(t)|y) -4(a;|r,(t)2G„(t)|?/) 
-4A J2 {^bix, tf + Gi,{x, y, t)) Gaix, y, t) . 



(3-6) 



Let us consider the case where the system has translational invariance. Then, the con- 
densates do not depend on space-coordinate x. Further, the two-point functions such as 
Ga{x,y,t) depend on the difference x — y only. Then, we can carry out the Fourier trans- 
formation in the same way as those in the section 2. Comparing Eq.(2-10) with Eq.(3-6), we 
here regard A<^^)[^] inEq.(2-10) as TW^^^ [^^j +4A I]fe^a[^6 + Performing the mode ex- 

pansion similar to Eqs.( |2-12| ) and ( |2-14| ), we finally obtain the following equations of motion 
for the condensate (mean field) and the quantum fluctuations : 



(2^ 



'{ty\if,{t)-h6, 



aO 



(2^ 



(3-7) 



These are the basic equations in the linear sigma model in the squeezed state approach. The 
condensate and the quantum fluctuations are coupled each other and the time-evolution of 
both the degrees of freedom should be determined self-consistently. 

It is instructive to deal with the static case in the 0(4) linear sigma model. The differen- 
tial terms with respect to time t disappear in this case. Let us assume that = = ^3 = 0^ 
that is, the chiral condensate points in the sigma-direction. Also, it is assumed that the each 
fluctuation mode of the direction of pi is identical, which we denote r]'^ : ~ '^k ~ ~ Vk- 
Then, the equations (3-7) read 



^Pq — h = 



+ Ml 



(3- 
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where we define the sigma meson mass M„ and the pion mass M^^ as 

Ml = -m^ + UXipl + 12A vi' + 12A r^£,^ , 

This is identical with the results derived in the previous paper. 0^ (The parameters A/24 
and ml should be replaced into A and — m^, respectively, in this paper. ) From Eq.(3-8), 
?7^^ and i]^^ are expressed as 

= -r== ' = / / (3-10) 



and the two point function Gq and are also obtained as 

^0 = [ /„ NO ; ^ ; G^TT = / /„ NO ; ^ • (^'H) 



The equations (3-8)~( |340| ) compose the set of self-consistent equations. 



§4. Numerical result and discussion 

We investigate the time-evolution of the mean field and quantum fiuctuations around 
it based on the above-derived equations of motion in the context of the quench scenario 
in the relativistic nucleus- nucleus collision as was mentioned in §1. Our attention is paid 
to describe the damping of the amplitude of the mean field. It will be shown numerically 
that the fluctuation modes are responsible for the damping of the mean field configuration. 
Namely, since the energy stored in the mean field configuration fiows to the fiuctuation 
modes, the behavior like the damped oscillation appears in the mean field configuration. 
Our aim is here to describe the damping behavior qualitatively. 

In general, the renormalization procedure is necessary in our approach because the loop 
contribution is included through the function G{x, x, t). The renormalization will be carried 
out because the linear sigma model is renormalizable. However, it should be noted that the 
sigma model is here regarded as a low energy effective model of QCD. Thus, we become free 
from the complexity of the renormalization due to introducing the momentum cutoff in the 
context of the low energy effective theory. 

The model parameters are given so as to reproduce the pion mass (M,r = 138 MeV), the 
sigma meson mass {M^ = 600 MeV) and the pion decay constant (/^ = 92 MeV) in the 
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condensate (signna,n=8) 




Time (fm/c) 

Fig. 1. The time-evolution of the chiral condensate is depicted in which the fluctuation modes are 
contained up to = 8^. 



static case, respectively, that is. 



Ml - Ml 



h 



8(/2 + Go - G^) ' 
+ Co — Gjr 



(4-1) 



where Gq and have been defined in Eq.( ^-llD and can be calculated with an appropriate 
three-momentum cutoff. In general, the initial condition of the mean field is taken arbitrary 
in our framework. However, the quantum fluctuation modes should be taken so that the 
quantum effects are as small as possible in static case. Thus, we adopt the initial conditions 
of the quantum fluctuation modes so as to satisfy dEa/dr]^ = with = 92 MeV, where 
represents the energy of the quantum fluctuation part. Also, we take rj^it = to) = 0. 
Let us investigate the time-evolution of the mean field on the a-direction only. We assume 
= with i = 1 ~ 3. In numerical calculation, we adopt the box normalization with the 
spatial length being L for each direction. We then impose the periodic boundary condition for 
the fluctuation modes, namely, the allowed values of momenta are = {2'k/ Vju^ and so on, 
where Ux is integer. The fluctuation modes labeled by (nj., Uy, Uz) are included in each isospin 



direction up to 



8^. This corresponds to the momentum cutoff /I ~ 1 GeV 



(990 MeV) since we have adopted the coUisional region as = 10^ fm^. It should be here 
noted that the equations of motion in Eq.(3-7) has the time-reversal invariance. Thus, the 
amplitude of the oscillation of the chiral condensate becomes large again even if the behavior 
of the damped oscillation is realized. However, in the realistic situation, the quantized pions 
and/or sigma mesons are emitted with a certain energy. Therefore, in our approach, it is 
enough to investigate the behavior of the condensate before the amplitude increases again. 
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Fig. 2. The time-evolution of the fluctuation modes in the cr-direction are depicted. The figures of 
the upper-left, upper-right, lower- left and lower-right represent the time-evolution of the fiuc- 
tuation modes with the quantum numbers = 0, = 2'^ , n? = 4^ and = 8^, respectively. 



Figure 1 shows us the time-evolution of the chiral condensate with Tp(t = to = 0) = 50 
MeV. It is shown that the behavior of the damped oscillation is realized. It should be noted 
that the center value of oscillation is slightly different from 92 MeV. The reason is that the 
effective potential contains the effects of quantum fluctuations in our approach. If 

the quantum fluctuation modes have no time-dependence, the value of condensate should be 
92 MeV. It is seen in Fig. 2 that the cr-modes oscillate. However, in Fig. 3, it is seen that the 
amplitude of 7r-modes with low momenta {n^ < 2^) increases. Namely, these amplitude of 
TT-modes with low momenta increases corresponding to the decreasing of the amplitude of 
the chiral condensate. However, the 7r-mode with high momenta are not amplified but oscil- 
lated only. This situation can be easily understood. The instability occurs when the time- 
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Fig. 3. The time-evolution of the fluctuation modes in the 7r-direction are depicted. The figures of 
the upper-left, upper-right, lower- left and lower-right represent the time-evolution of the fiuc- 
tuation modes with the quantum numbers = 0, = 2^, = 4^ and = 8^, respectively. 



dependent "mass" term Maiff' = — + 12A^^(t) + I]b^a 4A^f,(t) + (quantum correction) 
is negative and the absolute value of it is greater than k^. For a-mode, the time-dependent 
mass is Mo(t)^ = —m'^ + 12A^o(^) + (quantum correction) as is similar to Eq.(3-9). However, 
for the TT-mode, the mass M^^lt)'^ is less than Mo(t)^ because Mn-(t)^ = —m^ + 4A^Q(t) -|- 
(quantum correction). Thus, the instability occurs in the 7r-modes with low momenta. As a 
result, the behavior of the amplified oscillation is realized. On the other word, as was shown 



in §2, if the mode functions are expressed as ?7k(t) = ^/y'^^kit): the instability is realized 
when e^(i) < 0. When e^(i) < for |fc| < ■\/—M'^, then G{x,x,t) is obtained as 

G{x^x^t)= E ^^^^'^*'^'^'''''+ E (4-2) 
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Fig. 4. The time-evolution of the chiral condensate is depicted in which the fluctuation modes are 
contained up to = 4^ . 
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n 



The time-evolution of the fluctuation modes with the quantum numbers = (left) and 
= 3^ (right) are depicted in the cr-direction. 



Namely, this means that the imaginary part appears in the equation of motion for the mean 
field ^a(t). In this case the amplitudes of the low momentum modes with < yJ—M"^ of 
quantum fluctuation increase in time. Thus, the relaxation process for the chiral condensate 
occurs due to the instability of the quantum fluctuation. Then, the energy stored in the 
chiral condensate flows to the quantum fluctuations. 

In Figs.4~6 and 7~9, the only a few fluctuation modes up to — 4^ and = are taken 
into account in each isospin direction, respectively. The time-evolution of the condensate, the 
fluctuation modes with u-direction and the fluctuation modes with 7r-directions are depicted 
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Fig. 6. The time-evolution of the fluctuation modes with the quantum numbers = (left) and 
■nP = 3^ (right) are depicted in the 7r-direction. 
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Fig. 7. The time-evolution of the chiral condensate is depicted in which the fluctuation mode with 
n = is only contained. 



in Figs. 4, 5 and 6, respectively, up to = 4^. The qualitative behavior is almost same as 
that in Figs.l~3. The fluctuation modes with > 3^ for vr-mode are not responsible for the 
damping of the amplitude of the condensate. However, in Figs.7~9, the qualitative behavior 
is quite different from Figs.l~3 and 4~6. The condensate is relaxed to the vacuum value no 
longer. The reason why the relaxation process is not realized is that the quantum fluctuation 
modes are not included sufficiently. Thus we conclude that the quantum fluctuation mode 
for TT-direction with low momenta are responsible for the relaxation process of the chiral 
condensate during the chiral phase transition. 
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Fig. 8. The time-evolution of the fluctuation 
mode with n = in the cr-direction is de- 
picted. 



Fig. 9. The time-evolution of the fluctuation 
mode with n = in the 7r-direction is de- 
picted. 



condensate (sigma,n=8) 




Fig. 10. The time-evolution of the chiral condensate on the expanding geometry is depicted in 
which the fluctuation modes are contained up to = 8^. 

Let us consider the expanding geometryil) to break the time-reversal invariance. The 
basic equations are given in Appendix. The expanding parameter a{t) is assumed so as to 
be proportional to time : a{t) = t/ro + 1. The parameter is taken as tq = 1 fm/c. The other 
parameters are same as those in Figs.l~3. Also, it is assumed that the expansion occurs in 
D = 3-dimension. Figure 10 shows us that the behavior of the damped oscillation appears 
in the chiral condensate and the oscillation converges to the vacuum value. It is also seen 
that the amplitude of the quantum fluctuation modes also decreases gradually in time in 
Fig.ll (a-modes) and Fig. 12 (vr-modes). 



§5. Summary 

In this paper, we have investigated the time-evolution of the chiral condensate and the 
quantum fluctuation around it in the context of the chiral phase transition from the symmet- 
ric to symmetry-broken phase. We have formulated the time-dependent variational method 
with the squeezed state for 0(N) linear sigma model. It has been shown numerically that 
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The time-evolution of the fluctuation modes with the quantum number n 
4? (right) in the cr-direction on the expanding geometry are depicted. 



(left) and 



the behavior hke the damped oscillation for the chiral condensate appears due to the in- 
clusion of the quantum fluctuations in the squeezed state approach. We have pointed out 
that the quantum fluctuations, especially vr-modes with low momenta, are responsible for 
the occurrence of the relaxation process of the chiral condensate. Of course, the expansion 
of the system leads to a friction term.i'''0) However, even if the expansion is not taken into 
account, it is concluded that the relaxation process occurs due to the quantum fluctuations 
around the mean fleld conflguration. It may be also interesting to investigate the implication 
to the parametric resonance in the vr-mode. 



Acknowledgements 

One of the authors (Y.T.) would like to express his sincere thanks to Professor 
Dominique Vautherin for giving his interest in this work and stimulating discussions when 
he stayed in the Laboratoire de Physique Theorique des Particules Elementaires (LPTPE) 
of Universite Pierre et Marie Curie (Paris VI). He also wishes to thank Dr. H. Fujii in 
University of Tokyo for the early collaboration in this work and giving him useful comments. 



16 



fluctuation (pi,n=0/8) 




4 6 

Time (fm/c) 



fluctuation (pi,n=4/8) 



Fig 



. 12. 

r? - 




The time-evolution of the fluctuation modes with the quantum number n 
4? (right) in the vr-direction on the expanding geometry are depicted. 
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Appendix A 
Expanding geometry 



The collisional region may expand spatially. If it is necessary to include the effect of the 
expanding geometry, iil) the metric may be taken as 



ds'^ = dt^ - a{tfdx'^ = g^^{t)dx^dx'' 



(A-l) 



where a{t) represents the expansion of collisional region, which is governed by another dy- 
namics. Here, a{t) is given by hand. Hereafter, it is assumed that the expansion occurs in 
D = 1- oi D = 3-dimension. Further, in D = 1, the direction of the expansion is adopted as 
the z-direction. The metric tensor g^u{t), its inverse and its determinant g are given by 
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-l/a(t)V 



(A-2) 



with D = 1 or 3. Then, equations of motion in Eq.(3-7) are modified because the d'Alembertian 
d"^ /dt^ — V^, which acts on ^„ and is replaced by 



I d_ 

d_ 
di 



'-g 



-gg 



dt 



1 



-g^TT^^^a 



df" a{t) dt 
Namely, the d'Alembertian is replaced into 



(A-3) 



(A4) 



where the variable with suffix D (= 1,3) is defined as 

1 .„o „o. 1 



(f + f 

\D-\ ^•'x ^ Jy 



-f 

I J z 



(A-5) 



As a result, the equations of motion we solve are modified from Eq.(3-7) as follows : 



{tf Tp^it) - hSaO = , 



kl-m' + 12Xlp^(t)' + 12A y — r)^,(t) 



a u\2 



(2nf 



6^0 \ 



Sk' 
(2^ 



(A-6) 



Of course, the term with the first derivative of time causes the friction and the dissipative 
process is realized. 
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